Interplay of mesoscopic and Kondo effects for transmission amplitude of few-level 

quantum dots 
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The magnitude and phase of the transmission amphtude of a multi-level quantum dot is calculated 
for the mesoscopic regime of level spacing large compared to level width. The interplay between 
Kondo correlations and the influence by neighboring levels is discussed. As in the single-level case, 
the Kondo plateaus of magnitude and phase disappear with increasing temperature. At certain 
gate voltages, "stationary" points are found at which the transmission phase is independent of 
temperature. Depending on the mesoscopic parameters of the adjacent levels (like relative sign and 
magnitude of tunneling matrix elements), the stationary points are shifted to or repelled by the 
neighboring level. 
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I. INTRODUCTION 

In a remarkable series of experimentsji'^'^'^'^'^ the 
Heiblum group has analyzed the complex transmission 
amplitude, td = of a quantum dot embedded in 

an Aharonov-Bohm ring. In particular, by analyzing the 
Aharonov-Bohm oscillations of the conductance of such 
a ring, the dependence of both the magnitude and phase 
of the transmission amplitude, \td\ and a, were measured 
as a function of various parameters such as gate volt- 
age Vg applied to the dot, temperature T, mean coupling 
strength to the leads F, etc. 

The first two experiments in this series^i^ dealt with 
large dots containing many (> 100) electrons. The exper- 
iment by Yacoby et alX showed that coherent transport 
through a quantum dot is possible despite the presence of 
strong interactions. The next experiments by Schuster et 
al^ generated tremendous interest because the behavior 
of the transmission phase showed a surprisingly "univer- 
sal" behavior as function of gate voltage: the phase ex- 
perienced a series of sudden jumps by — tt (phase lapses) 
between each pair of Coulomb blockade peaks in the con- 
ductance through the dot. This contradicted a naive ex- 
pectation that the behavior of the transmission phase 
should depend on microscopic details of the dot, such as 
the signs of the matrix elements coupling a given level to 
the left or right lead. 

Subsequent experiments by Ji et al.^'^ performed on 
smaller dots containing tens of electrons, analyzed how 
the occurrence of the Kondo effect influences the trans- 
mission amplitude, and in particular its phase. For 
transmission at zero temperature through a single level, 
the Kondo effect causes the magnitude of the transmis- 
sion amplitude to exhibit (as function of gate voltage) a 
plateau at the unitary limit {\td\ = 1). For this regime 
it had been predicted by Gerland et al^ that the phase 
should show a plateau at a = 7r/2, a result very different 
from the universal behavior mentioned above. While the 



experiments of Ji et al. did yield deviations from the uni- 
versal phase behavior, they did not verify the prediction 
of a 7r/2 Kondo plateau in the phase. With hindsight, the 
reason probably was that the experiments did not realize 
the conditions assumed in the calculations of Gerland et 
al.^ namely transport through only a single level. 

Truly "mesoscopic" behavior for the phase was ob- 
served only rather recently by Avinun-Kalish et al.^ in 
even smaller dots containing only a small (< 10) num- 
ber of electrons. For these, the mean level spacing S 
was significantly larger than the average level width F, 
so that for any given gate voltage, transport through the 
dot is typically governed by the properties of only a sin- 
gle level, namely that closest to the Fermi energies of 
the leads. When the number of electrons was increased 
beyond about 14, universal behavior for the phase was 
recovered. Consequently, it was proposed^i^ '^^i^^'^^i^^ 
that the universal behavior occurs whenever a quantum 
dot is large enough for that the ratio 6/T is sufficiently 
small (~ 1) that for any given gate voltage, typically 
more than one level contributes to transport. 

The latest paper in this series, by Zaffalon et al.^ stud- 
ied the transmission phase through a quantum dot in the 
"deep mesoscopic" regime 5/T » 1, containing only one 
or two electrons. When this system was tuned into the 
Kondo regime, the transmission phase indeed did show 
the 7r/2 Kondo plateau predicted by Gerland et alX. 

The experiments of Avinun-Kalish et al.^ which ob- 
served mesoscopic effects for the transmission phase 
through a small number of levels, and those of Zaffalon et 
al.^ which found characteristic signatures of the Kondo 
effect in the transmission phase through a single level, 
raise the following question: what type of phase behav- 
ior can arise in the deep mesoscopic regime from the in- 
terplay of (i) random signs for tunneling amplitudes of 
neighboring levels and (ii) the Kondo effect for individual 
levels? In the present paper, we address this question by 
studying spin-degenerate models of dots with 2 or 3 lev- 
els in the deep mesoscopic regime of 5/F ^ 1. This is the 
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regime relevant for the experiments of Zaffalon et al^ (for 
those of Ji et al.^^ the ratio 5/T was presumably smaller 
than used here). Our goal is to provide a catalogue of the 
types of behavior that can occur in this regime, and to 
illustrate how the characteristic transmission amplitude 
(magnitude and phase) depends on temperate as well as 
on the strength of the coupling to the leads. 

This paper is organized as follows. In Sec. wc in- 
troduce our many-level model for the quantum dot sys- 
tem. We discuss the relation between the Aharonov- 
Bohm contribution to the linear conductance and the 
transmission amplitude through the quantum dot. The 
latter can be expressed in terms of the local Green's func- 
tion of the dot. We briefly present the technique used to 
calculate the latter, the numerical renormalization group 
method. In Sec. IIIII we present our numerical results of 
both the phase and the magnitude of the transmission 
amplitude through a two- and three-level model in the 
regime 5/T ^ 1. We discuss the T- and F-dependence of 
the transmission amplitude with focus on the influence 
on Kondo correlations. We study all relevant choices of 
the mcsoscopic parameters given by the relative signs of 
the tunneling amplitudes of adjacent levels. The influ- 
ence of neighboring levels is studied. It results not only 
in a phase lapse in Coulomb blockade valleys but also in- 
troduces a Vg-asymmetry in the finite temperature mod- 
ulations of the Kondo plateaus. "Stationary" points of 
T- and F-independcncc are discussed. In the Appendix, 
we give a derivation of a formula for the Aharonov-Bohm 
contribution to the linear conductance through a multi- 
terminal interferometer with open geometry, as used in 
the Heiblum group. This formula has been used in several 
publications including some of the present authors fi^iii^ 
but its derivation had not been published before. 



II. THE MODEL AND THE METHOD 

In the expcrimcnts^i^i^i^i^ the temperature-dependent 
transmission amplitude through the quantum dot is ex- 
tracted from the Aharonov-Bohm oscillations of the con- 
ductance in a multi-lead ring geometry. In the Appendix 
we show that this transmission amplitude can be ex- 
pressed in terms of the equilibrium local Green's func- 
tion of the dot tunnel-coupled only to two leads on its 
left and right side, without explicitly incorporating the 
other leads of the ring geometry in the calculation. 

In this Section we introduce a "reduced model" de- 
scribing the latter situation of a spinful multi-level quan- 
tum dot coupled to two reservoirs and present the trans- 
mission formula derived in the Appendix. Further, we 
comment on NRG, the method used to calculate the lo- 
cal Green's function. 



1. The model Hamiltonian 

The model Hamiltonian can be split into three parts, 

H = Hd + Hi + Ht, (la) 

specifying the properties of the bare dot, the leads and 
the coupling between the two systems, respectively. For 
N spinful levels coupled to a left (emitter) and right (col- 
lector) lead, these terms arc given by 

Hd ^ ^ ^£djndja+ ^ Uudjandj'a' (lb) 
(x=L,R ka 

= T.T. E(^'c.cL.rf..+H.c.). (Id) 

j a — L,R ka 

Dot creation operators for level j and spin (t = {|,|} 

are denoted by dj^, with Udja = dja^ja, where j = 
1 ■ ■ ■ N labels the levels in order of increasing energy 
{sdj < Edj+i)- We use an inter- and intra- level inde- 
pendent Coulomb energy U > 0. The leads are as- 
sumed to be identical and non-interacting with constant 
density of states p = 1/2D, where the half-bandwidth 
D = 1 serves as energy unit. Electrons in lead a are cre- 
ated by c|jj,^. The local levels are tunnel-coupled to the 
leads, with real overlap matrix elements that for sim- 
plicity we assume to be energy- and spin- independent. 
The resulting broadening of each level is given by Tj = 
TjL + Fjjj, with Tja ~ np{ti^)^. Notation: We define 
Si = sgn{t]^t]^t^j^^fj^^) ~ ±. For example, matrix ele- 
ments of same sign result in Si = +, whereas one different 
sign yields Si = —. We further define s = {si ■ ■ • sn-i}, 
and use 7 = {Fil, Fi_r, • • • , Fatl, Fat^rI/F, with the 
mean level broadening F ~ l/N'^j Tj. We assume con- 
stant level spacing S — Sdi+i ~ £di- The local levels can 
be shifted in energy by a plunger gate voltage Vg, with 
£dj = jS - {Vg + Vgo), where Vgo = ^6 + ^U. This 
convention ensures that in case of maximal symmetry 
(t^ = const, for all j,a), the system possesses particle- 
hole symmetry at Vg = 0. 

2. Transmission 

In the Appendix we generalize a result of Bruder, Fazio 
and Schoellei— to show that the Aharonov-Bohm con- 
tribution to the linear conductance through the multi- 
terminal interferometer with open geometry with a multi- 
level quantum dot embedded in one arm (see Fig. [3] in 
the Appendix) can be expressed as 

gAB(y) ^ El\T^\\t^{T)\ cos(2^$/$o + <^o + «(r)). (2) 
ri 

Here r„ = |r„|e*'^»+*2''*/*« is the energy- and 
temperature-independent transmission amplitude 
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through the upper reference arm including the Aharonov- 
Bohm contribution 27r<i>/$o to the phase, where $ is the 
magnetic flux enclosed by the interferometer arms and 
$0 = h/e is the flux quantum. The equilibrium Fermi 
function of the leads are denoted by /q. The effective, 
temperature-dependent transmission amplitude td{T) 
through the lower arm including the quantum dot is 
given by 

U{T) =JdE T,(i?,T) HE \t,\ e", (3) 

where 

UE,T) = ^^27rp titig^„^^,^,iE,T). (4) 

jj' o-ct' 

Therefore, only local properties like the local retarded 
Green's function Q^^^,^, and the Fermi function of the 
leads enter in the transmission amplitude through the 
quantum dot td [Eq. ([3])]. Thereby the local Green's 
function is evaluated for the model given in Eqs. ([T]) in 
equilibrium at temperature T. 

In the zero temperature limit and in linear re- 
sponse, the dot produces purely elastic potential scat- 
tering between left and right leads, which can be fully 
characterizedi^ by the eigenvalues e'^"'" {ly—a, b) per spin 

of the S'-matrix, and the transformation | ^ „ ^^'^ ^ I , 

smt^ coaO J 

that maps the left-right basis of lead operators onto the 
a-b eigenbasis of S. The transmission amplitude through 
the dot then reads 

td = -i5Lfl = sin(20) sin(<5, - ^,,)e*(''"+''''), (5) 

where in general 9 and Si, are all 1/g-dependent. The 
phase Si, is related by the Friedel sum rulei^ to the charge 
(per spin) rii, = Snj'K extracted by the dot from effec- 
tive lead V. As Vg is swept, the transmission amplitude 
goes through zero whenever = riftmodl, and a phase 
lapse by tt occurs. Equation ^ is useful for the spe- 
cial case of "proportional couplings", = i^^]?, with 
A independent of j, in which the occupations Ua.h take 
a simple form. Then the two effective leads a and h 
are the even and odd combinations of the left and right 
leads, respectively, with tan 6* = 1/A independent of Vg. 
Then each level either couples to the even or the odd 
lead, and the occupations extracted from the leads are 
given by UEfi — J2jeE o^dja-- Note that if all levels 
are coupled to the same effective lead (which is the case 
for s = {+■••+}), the other effective lead decouples, 
thereby reducing the computational complexity signifi- 
cantly. 

5. The method 

We calculate the local Green's function needed 
for the transmission amplitude (Eqs. ^ and re- 
spectively) using the numerical renormalization group 



methodic (NRG), a well-established method for the 
study of strongly correlated impurity systems. For a re- 
view, see Ref. [l^ The key idea of NRG is the logarithmic 
discretization of the conduction band with a discretiza- 
tion parameter A > 1. As a result. Hi is represented as 
a semi-infinite chain, where only the first site couples to 
the local level. The hopping matrix elements along the 
chain fall off exponentially like A~("~^)/^ with the site 
number n (energy scale separation). The NRG Hamil- 
tonian can be solved iteratively by successively adding 
sites and solving the enlarged system, thereby increasing 
the energy resolution with each added site by a factor 
of A^/^. The corresponding increase in Hilbert space is 
dealt with by a truncation strategy that keeps only the 
lowest iVkcop states for the next iteration. 

For the calculation of IniQ^ we use the full density 
matrix NR G^^i^° , based on the only recently developed 
concept of a complete basis set within NRG^. The real 
part of g^ is obtained by Kramers-Kronig transforma- 
tion. Improvement of the results is obtained by the self- 
energy representation, where the /7-dependent part of the 

impurity self-energy T,(uj) ~ U is expressed by two 

correlation functior^, which both are calculated with 
the full density matrix NRG. 

III. RESULTS 

In this Section we present our results for the phase and 
magnitude of the transmission amplitude td through the 
quantum dot. The gate voltage Vg is swept over a range 
sufficiently large that the full occupation spectrum of 
the quantum dot is covered ranging from to 2N. The 
exact distribution of the couplings seems to play only 
minor role for the transmission amplitude. Therefore 
we choose left-right symmetric coupling in the cases 
where all Sj = +, reducing the computational effort 
significantly, since then the odd channel decouples. 
In the regime of interest, the deep mesoscopic regime, the 
mean level spacing S is much larger than the typical level 
widths Fj, S/r ^ 1. Therefore electrons enter the dot 
one by one when increasing the gate voltage. Transport 
thus occurs mainly through one level at a time; more 
precisely, it occurs through a linear combination of all 
levels, where in the mesoscopic regime the level closest 
to the Fermi energy dominatesi^ 

The Section is organized as follows: We first cludicate 
the basic properties of the transmission amplitude for the 
example of a two-level system. Varying temperature T 
(at fixed coupling F), or average coupling F (at fixed T), 
we study both possible choices s = + and s = — , respec- 
tively. In order to analyze the interplay of s = -I- and 
s = — , we then present data for a three- level system for 
all four possible combinations of si, S2- Additionally, this 
has the advantage that for the middle level "boundary ef- 
fects" (effecting the outermost levels) can be assumed to 
be eliminated, thus the behavior of the middle level can 
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N=2, T-dependence 8/U=0.3, yg^^H0.4,0.4,0.6,0.6}, y^^ ={0.3,0.5,0.7,0.5}, U=0.6,A=3.2 




FIG. 1: Transmission td = ltdle'™ through a spinful two-level quantum dot for various temperatures and constant couplings. 
Regimes (i), (ii), indicated in panel (a) only, refer to Kondo valleys or Coulomb blockade valleys, respectively (see text). The 
levels involved are indicated by their level number 1,2. Level 2 is coupled more strongly to the leads than level 1, resulting 



in different bare Kondo temperatures, e.g. 



We use r/C/ = 0.03 (a,c) and T/U = 0.08 (b,d). 



thus r^'''''^< 



(b,d) 



The minimum value of the (in the center of the Kondo valleys) are indicated in the legends. In accordance with Ref. |23| . 
we find shoulders in the phase (see e.g. the fat purple arrow and the purple curve {T/U = 10~^) in (c) for level 1) and an 
enhanced sensitivity of the phase to Kondo correlations compared to the magnitude, see e.g. the green curves m 

(d) for level 1 or the purple curve {T/U = 10"*^) for level 2 in (a). There, the typical ^-Kondo plateau in the phase is present, 
whereas the Kondo plateau in amplitude is not fully developed yet. At certain points in gate voltage, say Vg^ (as indicated 
by red arrows), we find stationary points where the curves for a for all temperatures intersect. The position of Vg^ is shifted 
by the presence of a neighboring level, being repelled by or shifted towards the latter for s = + or — , compare (a,c) or (b,d), 
respectively. Depending on the mesoscopic parameter s = ±, the phase either exhibits a sharp drop of n, accompanied by a 
zero in the amplitude \td\ {s = +, see (a,b)), or increases monotonically (s = — , see (c,d)) in the Coulomb blockade valleys. 



be viewed as representative of a generic level in a multi- 
level quantum dot in the deep mesoscopic regime. 

Unless otherwise noted, we use U = 0.6. In order to 
cover all relevant energy scales with reasonable compu- 
tational effort, we usually use A = 3.2 for the two-level 
model and A = 3.5 in case of three levels. We checked 
that already by keeping ^ 1000 states at each iteration, 
also for the two-channel calculations (that involve at least 
one Si — ~) the physical trends are captured qualita- 
tively. Note that since the eigenvalues of the scattering 
matrix are given by e*^''" , the transmission phase a is de- 
fined modulo TT. For clarity of the Figures, curves showing 
a are shifted by multiples of tt as convenient. 



A. Two-level model 

1. Temperature dependence 

Figure [T] shows the transmission amplitude for both 
s = + (a,b) and s = — (c,d), for fixed dot parame- 
ters and various different temperatures. The mean level 
broadening is chosen to be F/C/ = 0.03 in panels (a,c), 
and T/U = 0.08 in panels (b,d). Therefore the {Vg- 



dependent) bare Kondo temperatures 

^^ = v^'"pr ^^^J (6) 

vary in a lower-lying range of energies for panels (a,c) 
than for panels (b,d). In all panels the relative couphng 
of the first and the second levels are chosen to be 
7 — {0.8, 1.2}. Therefore, the bare Kondo temperature 
for level 1 is lower than for level 2, < as 

indicated in the legends. The resulting difference in the 
temperature dependence can be nicely observed in the 
Figure. We first describe those general properties of the 
transmission amplitude that qualitatively agree with 
those that one would obtain for just a single level, then 
discuss the effect of the presence of a second level. 

General properties: In the mesoscopic regime, where 
transport mainly occurs through one level at a time, two 
different regimes of transmission can be distinguished as 
Vg is varied, as indicated in Fig. [IJa): (ii) In the regime 
between the Kondo valleys, to be called "Coulomb block- 
ade valleys" , the transmission amplitude is mainly deter- 
mined by the mesoscopic parameter s, showing a phase 
lapse only in case s = +, similar for both spinful and 
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spinless model a^^i^^ . 

In the zero temperature limit, T <^ Tj^ , the transmis- 
sion amplitude exhibits the typical Kondo behavior: in 
the local- moment regime a typical Kondo plateau forms, 
with \td\ approaching the unitary limit, ^ 1. In 
the mixed valence regime the magnitude changes rapidly 
as a function of V^. In the Coulomb blockade valleys, 
transmission is suppressed by Coulomb interaction. The 
transmission phase increases by ~ 7r/2 for each entering 
electron (see black curves for a in Fig. [T|) , increasing only 
slightly in between. In the Kondo valleys this results in 
a plateau at amodTr = ^ , as direct consequence of the ^ 
phase shift due to the formation of the Kondo singlet. 

With increasing temperature, the Kondo effect is sup- 
pressed, thus the behavior in the middle of the Kondo 
valleys changes dramatically. The Kondo plateaus in td 
and a disappear: The magnitude tends towards Coulomb 
blockade behavior, with a resonance of width ~ V j for 
each entering electron. The phase develops a S-like shape 
in the Kondo valleys with increasing temperature. As in 
the single-level case, all finite-temperature curves of the 
phase intersect the zero-temperature at the same gate 
voltage, say Vg' (see red arrows). We shall refer to this 
gate voltage as a "stationary" point (w.r.t. temperature). 

As observed in the experiments of Ji et al^ and 
emphasized by Silvcstrov and Imry,— the transmission 
phase reacts more sensitively to the buildup of Kondo 
correlations with decreasing temperature than the trans- 
mission magnitude: a approaches its T = behavior 
already at temperatures T ~ Tk (the ^-plateau devel- 
ops), whereas \td\ develops its plateau for T significantly 
less than Tk (see, the green curve {T/U = 10~^) for level 

1 in Fig.[Hb) or the purple curve {T/U = 10~^) for level 

2 in Fig. [I^a)). Similar to the predictions of Silvestrov 
and Imry,— we find shoulders in the evolution of the 
phase, see for example the fat purple arrow and the 
purple curve [T/U = 10^^) in Fig. [T]Jc). This indicates 
that the temperature is large enough to suppress Kondo 
correlations in the deep local-moment regime (in the 
middle of the Kondo valley), where Tk is very small. 
Towards the borders of the local-moment regime the 
crossover temperature for the onset of phase sensitivity 
increases (as does the Kondo temperature, see Eq. 
([6])), eventually exceeding the temperature. Then the 
phase tends towards its zero-temperature behavior, thus 
producing shoulders. 

Properties special to the multi-level model: The most 
obvious difference between the transmission amplitude of 
the many-level model in the mesoscopic regime compared 
to the single-level model is the phase behavior in the 
Coulomb blockade valleys between the levels. Depending 
on s, i.e. on the relative sign of the tunnelling matrix ele- 
ments of the two adjacent levels, the phase either exhibits 
a sharp drop (phase lapse) by tt in the s = + case (ac- 
companied by a transmission zero, \tfi\ = 0), or evolves 
continuously for s ~ 11, 12, 14, 24, 25, 26 _ Contrary to the 
non-monotonic phase evolution discussed above, this ef- 



fect occurs already at zero temperature and also exists 
for spinless models i^^i^^ Therefore, the relevant energy 
scale for the temperature dependence of this phase lapse 
is not related to the Kondo temperature but to the level 
distance and width of the effective transport levels 1^ It is 
therefore not a relevant energy scale in the temperature 
range studied in this work. 

A further peculiarity for models with more than 
one level is the asymmetry (w.r.t. the center of the 
Kondo valleys) of the transmission amplitude in the 
local-moment regime at finite temperature, introduced 
by the mixing of neighboring levels. The asymmetry 
in phase can be characterized by the position of the 
stationary points, Vg' (indicted by red arrows in Fig. 
[1]). In case s = +, these points are repelled by the 
neighboring level, whereas they are shifted to the latter 
for s = — , compare for example Fig. [TJa) and (c) or (b) 
and (d). For ri/r2 ^ 1, the repulsion and attraction 
is enhanced or reduced compared to Fi = T2 for the 
level that is coupled less or more strongly to the leads, 
respectively. Clearly, in the limit of one decoupled 
level (effective one-level system), the stationary point 
of the other level is symmetric w.r.t. the corresponding 
Kondo plateau. The dips that form in the plateaus of 
the amplitude with increasing temperature develop a 
distinct asymmetry only for T 3> , for which they 
tend to shift towards the corresponding Vg' . This is 
consistent with the fact that as the phase drop in the 
Kondo valley gets sharper with increasing temperature 
and approaches a quasi-phasc lapse, the magnitude 
experiences a minimum, as for every complex function. 
Interestingly, the asymmetry in phase is the same for 
all temperatures, thus already at temperature T < Tk 
the phase "knows" in which direction (of Vg) the dip in 
magnitude will shift at higher temperatures. 



B. Dependence on the coupling strength 

In experiments, it is more convenient (and easier to 
control) to change the coupling strength between the 
quantum dot and the reservoirs than the temperature. 
Accordingly, Fig. [2] presents the transmission amplitude 
for various values of F, keeping the temperature constant. 
With decreasing F, the decrease of Tk together with the 
suppression of Kondo correlations is nicely illustrated. 
At fixed temperature T > T/^, the S-like shape of the 
phase evolution gets more pronounced and sharper with 
decreasing F. 

In the single-level problem, in addition to stationary 
points w.r.t. temperature, we also find stationary points 
w.r.t. F for td^ i.e. for magnitude and phase of the trans- 
mission amplitude. These occur at the outer flanks of the 
Kondo plateaus. Varying the mean coupling strength F 
at fixed 7, 5 and T in the two-level model, as shown in 
Fig. [21 these points can still be recognized (indicated by 
green arrows in (d)), even though the F-indepcndence 
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N=2, r-dependence y3,^={0.4,o.4,o.6,o.6}, y^^={o.3,o.5,o.7,o.5}, u=o.6, a=3.2 




FIG. 2: Transmission through a spinful two-level quantum dot for both choices of s = ± and various values of mean couplings 
r at fixed temperature T, level spacing 5 and relative couplings 7. Due to the mixing of the levels, no stationary points w.r.t. 
r exist, F-independence exist, see text and the green arrows in (d). 



is not perfect (within our numerical accuracy). We ex- 
pect that due to the mixing of the levels, also the level 
distance 5 has to be taken into account to recover these 
stationary points. Between the levels, near Vg/U « 0, 
another stationary point seems to occur. 



C. Three- level model 

Naturally, the question arises about the effects of sev- 
eral levels, with different choices of = ±, which is 
present only for models with more than two levels. As- 
suming that in the mesoscopic regime only neighbor- 
ing levels mix significantly, i.e. simultaneously influence 
transport, any local level of a quantum dot (except the 
lowest or highest one) can be represented adequately by 
the middle level of a three-level model. 

In Fig. [3] we present numerical data of a three-level 
model for all four possible combinations of s = si,S2 
and various temperatures. The second level is influenced 
by the effect of both si and S2, resulting in an effective 
enhancement or compensation of the asymmetry of the 
stationary point Vg^ of level 2, as discussed in Section 
nil A 11 Also the relative strength of the level-couplings 
(given by 7) has to be considered. In Fig. El^a), both s 
and 7 symmetrize the transmission curves of the middle 
level, whereas in panel [2]^b) 7 shifts Vg^ to positive Vg. 
In panels [3jc) and [S^d) both s and 7 tend to increase 
the asymmetry. 

Therefore, the transmission phase through a spinful 
quantum dot with Kondo correlations present has S'-like 
shape in the local-moment regimes at T ^ Tk- Analo- 
gously to experiments, we find an asymmetry of this S'- 



like shape. It is determined by both the relative strength 
7 and the sign s of the level couplings. 



IV. CONCLUSION 

In this paper we present temperature-dependent NRG 
calculations of the magnitude and phase of the transmis- 
sion amplitude through a multi-level quantum dot in the 
regime 5/T <C 1. Clearly, the Kondo correlations are 
suppressed with increasing temperature. The presence 
of neighboring levels results in a Vg-asymmetry in the fi- 
nite temperature modulation of the Kondo valleys. The 
asymmetry depends on the relative signs of the tunnel- 
ing matrix elements as well as on the relative couplings 
of the adjacent levels. Further, sharp phase lapses may 
occur between the levels. Studying a three-level model, 
the middle level can be understood as a representative of 
a generic level in a multi- level quantum dot. 

Throughout the paper, we deliberately focussed only 
on the deep mesoscopic regime, for which the results can 
be understood rather straightforwardly. The crossover 
into the regime 5/T ~ 1, which is certainly of interest 
too in order to understand the fate of Kondo physics 
in the universal regime, and which we believe to be the 
regime relevant for the experiments of Ji et al.^^ will be 
left as a subject for future studies. 
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N=3, T-dependence 



r/U=0.03, T|^{Y=0.8)/U=9 10"^ T^(1)/U=3 10'^ 1,^(1 .2)/U=2 10"'^ 




FIG. 3: Transmission through a three-level model for all four possible combinations of s = Si, S2 and various temperatures, for 
fixed r/f/ = 0.03. The two-channel calculations for (b,c,d) qualitatively capture the physical trends. The asymmetry in the 
Kondo valleys is determined by both s and 7. For convenience the figure legends for 7 display only the total relative coupling of 
each level. The minimal bare Kondo temperatures are indicated, (a) s — ++: 7 = {0.6, 0.6, 0.4, 0.4, 0.5, 0.5}. The case T/U = 

10"® is included only for this panel, (b) s = : 7 = {0.5,0.3,0.4,0.6,0.5,0.7}. (c) s = -I--: 7 = {0.4,0.8,0.4,0.4,0.3,0.7}. 

(d) s = -+: 7 = {0.5, 0.3, 0.3, 0.7, 0.7, 0.5}. 
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2]) is chosen real, where t^^ = {xi\j(j) (blue) is the ampli- 
tude to get from dot state \ja) of level j and spin a to 
point Xi on side i = L,Roi the dot, and A\^^ = {£aa\xi) 
(red) is the amplitude to get from point Xi to lead state 
\£aa) in lead a with energy e and spin cr, sec Fig. 2] 



APPENDIX A: CONDUCTANCE FORMULA 
FOR MULTI-TERMINAL GEOMETRY 

1. General case 

We generalize the current formula derived in by 
Brudcr, Fazio and SchocUer— for a single-level quantum 
dot embedded into one arm of an Aharonov-Bohm 
interferometer with two-terminal geometry to a multi- 
terminal geometry with a multi-level dot (as used in the 
Heiblum group^i^iii^>^) . 

Consider a iV-level quantum dot described by [Eq. 
(jlbp ] embedded in one arm of an Aharonov-Bohm inter- 
ferometer connected to M leads, as depicted in Fig. U) 
Each lead, and each arm connecting them, is assumed 
to support only a mode. The tunnelling between the lo- 
cal levels j = 1 • ■ • iV on the quantum dot and the leads 
a = 1 • • • M is described by 



.t 



B.C. 



(Al) 




Here t|„^ 



T.^=L,R^laKaa (mdicatcd in grccu in Fig. 



FIG. 4: Geometry of the multi-terminal Aharonov-Bohm in- 
terferometer with a multi-level quantum dot embedded in the 
lower arm. The different tunnelling amplitudes used in the 
text are indicated. $ is the magnetic fiux penetrating the 
interferometer. 

Following Biittikerj^ the current operator in reservoir 
a is given by 

ee' oa' 

(A2) 

where p, the density of states, is assumed to be constant 
and equal for each reservoir. The first term inside the 
bracket stands for the incident, the second term for the 
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reflected current in reservoir a, thus bi^aa = X]/3 '^aff'^sfSa, 
with the scattering ampHtude to get from lead P to 
lead a with energy e. Defining the lesser, retarded and 
advanced correlation functions 



-{al,{t')a,{t)) 



dE 

2TTh 



Ql,,{El (A3) 



^ -'-6{±{t-t')){[al,{t'),a,{t)\^) 

- l^e-^^''-''/^'Of;;^{E), (A4) 

where denotes a fermionic operator with composite 
index fi, the expectation value of the current operator 
[2|) can be expressed as 



ee' PIS' aa' ^ 



(A5) 



To calculate Q^{E) in Eq. (|A5|) . we use the standard 
Dyson equation for the Keldysh 2x2 matrix Green's 
fimction2^ g{E), 

{E) = 5,e'5a.c'5a.'Q'',^AE) (A6) 

which yields 

Qfo..,e'.'AE) = 5,,,5o.^,5„„,g^<,{E) (A7) 
n' 

where terms in square brackets are given by 

A = g°^AE) Ql^y..{E) Q%AE) , 

B = g',t{E) g<,.,.{E) g%AE) , 
c = g'<AE) g^^,,,,{E) g%AE) , 

where the free Green's functions for the leads have the 
form 



(A8) 



gi<,{E) = 27:zUE)S{e ~ E), (A9) 

with fa{E) the Fermi function of lead a. Inserting g^ 
[Eq. (|A7p ] into Eq. (|A5p . the current can be written as 



(/„) = /::+5/„. 



(AlO) 



/° arises from the first term of Eq. (jA7p . It describes 
the situation when the dot is completely decoupled 
(^ectcr = 0); thus does not contribute to Aharonov-Bohm 



oscillations. The influence of the quantum dot on the 
Aharonov-Bohm oscillations is caused by (5/a, arising 
from the second expression of Eq. (|A7[) . Performing the 
energy sums X^ee' '^^ Eq. (jAlOp or (|A5|) , respectively, the 
two contributions to the current read 

= J / rfi?EE['5"/3-|^<^/3l']/M^)' (Alia) 
61^ = ^Rej /rfi?EEE['^"/3'<^a/3-5:f.5f;3] 
>^^PiW%'A-^) (Allb) 



^[2gf^,y^^{E)fp,{E) + g<^^,AE)] 



2. Simplification to effective 2-lead geometry 

For the experimental setup used by Schuster et al^ 
(and equivalently for the ensuing papers^i^i^i^) to mea- 
sure transmission phase shifts, two simplifying assump- 
tions can be made. The first allows us to neglect non 
equilibrium effects, the second to perform NRG calcula- 
tions for a simplified geometry, in which the dot is cou- 
pled only to two leads. 

(i) Neglect of non equilibrium effects: In the experi- 
mental setup used by Schuster et al.^ the leads a = 3, 4, 5 
and 6 serve as draining reservoirs (to prevent multiple 
traversals of the ring, see below), and are all kept at 
chemical potential = 0. This also fixes the chemical 
potential of the ring, referred to as "base region" in Ref. 

to equal zero. Lead 1 and 2 serve as emitter and col- 
lector, respectively, with chemical potentials fix and /j,2, 
and Fermi functions fi,2{E) = fo{E — A*i,2)- Now, the 
point contacts between emitter or collector and the base 
region (marked Pi and P2 in Fig. |4|) are so small that 
the voltage drops occur directly at these point contacts, 
and not at the tunnel barriers coupling the dot to the 
ring. Thus, while the emitter or collector inject or ex- 
tract electrons into or from the base region, respectively, 
this is assumed to happen at a sufficiently small rate that 
the base region is not disturbed. In other words, we may 
assume that the dot, ring, and electrodes 3, 4, 5, 6 are all 
in equilibrium with each other, and that the dot Green's 
functions g^^^,'^, {E) do not depend on jii and fj,2 at all. 
Thus, the lesser function can be expressed in terms of the 
retarded and advanced ones using the following standard 
equilibrium relation: 



g<,^.,,,(i?) = -fo{E) [gf,,,,,,{E) - gf,,,,,,{E)] . 

(A12) 

The conductance in the linear response regime can be 
obtained by taking fJ-i — fJ.2 ~ eV, where e = |e|, e.g. by 
setting 

Ml =0, = -eV, (A13) 

and calculating G = dli/dV, with Ii given by Eq. (jAlip . 



9 



(a) Reduction to two-lead geometry: The reason why 
a multi-lead geometry was used in experiment is to avoid 
phase-rigidity: in an Aharonov-Bohm ring connected to 
only two leads, the transmission phase of the dot does 
not vary smoothly with gate voltage, but can assume 
only two distinct values, differing by tt. A multi-lead ge- 
ometry avoids this by strongly reducing the probability 
amplitudes for paths from emitter to collector to traverse 
the ring multiple times, since with each traversal of the 
ring the probability increases that electrons travelling in 
the ring are "siphoned off" into the side arms. We shall 
exploit this fact by making the assumption that the prob- 
ability amplitudes for multiple traversals of the ring are 
negligibly small. This assumption allows us to replace the 
multi-lead geometry with one where the Aharonov-Bohm 
ring is coupled to only two leads, i.e. a is restricted to the 
values 1 and 2 (corresponding to emitter and collector), 
while multiple traversals of the ring are eliminated (by 
hand) by the following specification: The amplitude f^^^. 
to get from state |jcr) on the dot to state \eaa) in lead 
a is taken to be nonzero only for the short, direct path 
from the dot to lead a, without traversing the upper arm 
(more correctly: we take A^^^ = for a = 2,4, 6 and 
^eaa = ^'^^ ^ = IjSjS). When calculating the current 
we do allow for direct paths from lead 1 to 2 via the upper 
arm, and lump all flux-dependence into the correspond- 
ing scattering amplitude, taking g*^'^*/*". How- 
ever, the upper arm is ignored for the calculation of the 
equilibrium local retarded or advanced Green functions 
Qf^^i ^i{E) using NRG. For the latter purpose, we thus 
use a model of a multi-level dot coupled to two indepen- 
dent leads, say L and R, with equal chemical potentials 
/^L = M-Ri representing the two segments of the ring to 
the left and right of the ring, coupled to it by tunnelling 
contacts. (These two segments should be treated as in- 
dependent leads, due to the assumption of no multiple 
traversals made above.) With the assumptions (i) and 



(ii) just described, let us now obtain an expression for 
that part of the conductance showing Aharonov-Bohm 

oscillations with applied flux, G^^ — , where is 
that part of the current in lead 1 depending on e'^'^*/*". 
For the chemical potentials given by Eq. (|A13[) , this cor- 
responds to evaluate Eq. (|Allbp with a = 1, /3' = 2 and 
/3 = 1, and we readily obtain 

G^'^lr) ^^JdERe [T:iE)UE)] (-^^) , 

(A14) 

where 

w) = Y.T.'^^p^EiArAErL, (ais) 

T*{E) = z5i*f^fi = |T„(S)|e^(2^*/*"+'^"(^» (A16) 

may be interpreted as the transmission amplitudes 
through the lower and upper arms, respectively. 

Assuming the transmission amplitude T.^ through the 
upper arm to be energy- and temperature- independent, 
the Aharonov-Bohm contribution to the conductance is 
given by 

g2 

G^'^iT) ^ -^T„||trf(r)|cos(27r$/$o+0o+a(r)).(A17) 

Then, the temperature-dependent magnitude and phase 
of the transmission amplitude through the quantum dot, 

t,{T) = jdE {- ^^'ff^ ) UE,T) ^ |i,(T)|e-(^), 

(A18) 

can be (i) extracted via Eq. (|A17|) from the experimen- 
tal results as well as (ii) calculated with NRG using Eq. 
(IM) . 
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